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Abstract The X-ray dim isolated neutron stars (XDlNSs) are peculiar pulsar-like objects, 
r^ \ characterized by their very well Planck-like spectrum. In studying their spectral energy 

f^ I' distributions, the optical/UV excess is a long standing problem. Recently, Kaplan et al. 

(2011) have measured the optical/UV excess for all seven sources, which is understand- 
Jl^ . able in the resonant cyclotron scattering (RCS) model previously addressed. The RCS 

"^ ' model calculations show that the RCS process can account for the observed optical/UV 

Ci . excess for most sources . The flat spectrum of RX J2143. 0+0654 may due to contribution 

from bremsstrahlung emission of the electron system in addition to the RCS process. 
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|>! . 1 INTRODUCTION 

o. 

The seven X-ray dim isolated neutron stars (XDINSs) are puzzling pulsar-like compact stars. They are 
characterized by their very well Planck-like spectrum, and then they are good specimen for neutron 
star cooling, atmospheric and equation of state studies (Tong et al. 2007,2008; review TuroUa 2009). 
Additionally, these sources may also have some relations with other classes of pulsar-like objects, e.g., 
/\ ' magnetars and rotating radio transients, etc (Tong et al. 2010). 

The spectral energy distributions of XDINSs show that their optical and UV emissions are above the 
extrapolation of the X-ray blackbodies, i.e., optical/UV excess. This excess of XDINSs is a long standing 
problem which may be related to the equation of state of neutron stars (Xu 2002, 2009). Previously, only 
two sources (RX J 1856.5-3754 and RX J0720.4-3125) have both optical and UV data (van Kerkwijk & 
Kaplan 2007). Recently, Kaplan et al. (201 1) have presented results of the optical and UV emission for 
all seven sources. They also find that the source RX J2143. 0+0654 has a very flat spectrm {Fi, ex f", 
a ^ 0.5). This new observation requires updated modeling of existing theoretical models. 

Tong et al. (2010) have considered the resonant cyclotron scattering (RCS) in pulsar magneto- 
spheres as the origin of optical/UV excess in XDINSs, i.e., the RCS model. This model can explain 
the observations of RX J1856.5-3754 and RX J0720.4-3125 (Tong et al. 2010). Considering the recent 
observational progress, we apply the RCS model to all seven XDINSs in this paper 

Section 2 is summary of the RCS model. Its application to XDINSs is treated in Section 3. 
Discussions and conclusions are presented in section 4. 
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2 THE RCS MODEL: A SUMMARY 



Tong et al. (2010) have considered the RCS model for the origin of optical/UV excess in XDINSs. It 
has the following three key points. 



1. 



In analogy with magnetospheric study of magnetars and rotating radio transients, also similar to the 
"electron blanket" of neutron stars (Ruderman 2003), there may also be a dense electron system 
in closed field line regions of XDINSs. We call it the "pulsar inner radiation belt" (since it is not 
far away from the central star), in comparison with the "pulsar radiation belt" proposed by Luo & 
Melrose (2007) (It may be called the "pulsar outer radiation belt", since it is near the light cylinder 
far away from the central star). 

The scatterings between surface X-ray photons and the electron system are modeled three dimen- 
sionally using the Kompaneets equation method. Numerical calculations show that there exist both 
down scattering and up scattering of the RCS process. The existence of down scattering of RCS 
process is crucial for its application to XDINSs. 

The application of down scattering of RCS process to XDINSs can explain its optical/UV excess 
problem. 



The Kompaneets equation (photon diffusion equation) for RCS process is (Tong et al. 2010): 
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Here, n is the photon occupation number, its evolution with time (which corresponds to spectra change) 
is determined by the right hand side of equation ([T]i. .t is the dimensionless photon frequency in units of 
electron temperature, trcs is the RCS optical depth. In order to calculate the RCS modified blackbody 
spectrum, four input parameters are needed; the neutron star surface temperature, the temperature and 
density of the electron system, and a normalization constant. 

For photons in a certain frequency range, e.g., from optical to soft X-ray, only electrons in a specific 
radius range will contribute to the RCS process. The electron system may extend to much larger radius, 
however the electrons there will not contribute to the RCS process (Tong et al. 2010). On the other 
hand, the electrons will radiate thermal bremsstrahlung emission, which indeed depends on the radial 
extension of the electron system. The inclusion of thermal bremsstrahlung process in addition to the 
RCS process will result in a flat optical spectrum. 



3 APPLICATION OF THE RCS MODEL TO X-RAY DIM ISOLATED NEUTRON STARS 

The application of RCS model to RX J 1856.5-3754 is shown in figure[T] Model calculations for the rest 
six sources are shown in figure |2] The input parameters of RCS model are given is table 1 . The neutron 
star surface temperature T^ is chosen as the temperature of the X-ray blackbody (from table 2 in Kaplan 
et al. 201 1) and is kept fixed during the fitting process. The normalization is the radius of the radiation 
region over the source distance. For example, for RX Jl 856.5-3754, R/d = 10/0.16(km/kpc). It means 
that for a source distance 0.16 kpc the radius of the radiation region is 10 km. All distance data are from 
table 4 in Kaplan & van Kerkwijk (2009). The electron temperature To is about 0.5 T^ or lower. The 
electron number density A'c is around 1.5 x 10^^ cm^"^. The neutral hydrogen column density A^H,rcs is 
one and a half or two times that of blackbody fit to the X-ray data. Considering that for all seven XDINSs, 
their surface magnetic field, surface temperature, radiation region. X-ray luminosity, and distance are all 
similar to each other (Kaplan & van Kerkwijk 2009). Meanwhile, they all have an optical/UV excess of 
order ten etc (Kaplan et al. 201 1). Hence, these seven sources are similar to each other observationally. 
As a consequence of this, the RCS model parameters for these sources are also similar to each other 

The two sources RX J1605.3H-3429 and RX J2143.0H-0654 have flat optical/UV spectra (especiafly 
for RX J2143.0H-0654). For most sources, a power law fit to the optical/UV data gives Fi, oc z^", with 
a ^ 2 implying a thermal origin. While for RX J2143.0H-0654, a ^ 0.5 (table 5 in Kaplan et al. 
2011). The flat spectrum of RX J2143.0H-0654 indicates that nonthermal processes are involved in its 
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Fig. 1 Spectral energy distributions of RX J 1856.5-3754 and RCS model calculations. The 
gray region is the X-ray blackbody and its extrapolation to the optical and UV range, including 
1 -G uncertainties. The red region is the power law fit to the optical and UV data, including 1 -a 
uncertainties. The red circles are observational data points, only central values are included. 
All the observational data are from Kaplan et al. (2011). The gap between optical/UV and 
X-ray is due to neutral hydrogen absorption. The dashed line is RCS model calculation (Tong 
et al. 2010). 



optical/UV emission. Noting that the bremsstrahlung emission has a flat spectrum, we propose that the 
flat spectrum of RX J2143.0H-0654 may due to thermal bremsstrahlung emission of the electron system 
in addition to the RCS process. 

From table 1, the temperature and density of the electron system are known. Once its radial exten- 
sion is given, we can calculate its bremsstrahlung flux at earth (Rybicki & Lightman 1979). The outer 
radius of the electron system rout determines the relative importance of the bremsstrahlung component. 
For RX J1605.3H-3429 and RX J2143.0H-0654, rout is determined in order to obtain a flat optical spec- 
trum. For the rest five sources, no bremsstrahlung component is needed. Therefore, an upper bound on 
''out is given. The rotation periods of XDINSs are about 10 sec. Therefore, the radial extension of the 
electron system is at most five percent of the light cylinder radius, i.e., near the central star This is why 
we call it "pulsar inner radiation belt". 

When the electron number density varies with radius, e.g, Ne{r) ~ N^ x (r/_R,is)~", the fitting 
results will be different from the uniform density case as considered above. In order for electrons in 
a large space range all contributing to the bremsstrahlung process, the power index a should not be 
significantly larger than 1, i.e., < a < 1 is required for RX J1605.3H-3249 and RX J2143.0H-0654. 
While for the rest five sources, we do not see such a bremsstrahlung component. This is due to either a 
small radial extension of the electron system or the electron number density there decreases with radius 
very rapidly, i.e., a > 1. The formation of the electron system is discussed in the appendix. 
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Fig. 2 Calculations for RX J0420.0-5022, RX J0720.4-3125, RX J0806.4-4123, RX 
J1308.6+2127, RX J1605. 3+3249 and RX J2143.0+0654. The same as figure 2. All obser- 
vational data are from Kaplan et al. (2011). The dashed line is the RCS model calculation 
(Tong et al. 2010). The dotted line is the contribution from electron thermal bremsstrahlung. 
The dot-dashed line is the sum of bremsstrahlung and RCS component. The bremsstrahlung 
component is needed only in the case of RX J1605. 3+3249 and RX J2143.0+0654. 
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Table 1 Input Parameters of the RCS Model. Column one to seven are source name, neu- 
tron star surface temperature, normalization, electron temperature, electron density, neutral 
hydrogen column density, and outer radius of the electron system, respectively. See text for 
details. 



Source 


(eV) 


R/d 
(km/kpc) 


(eV) 


(10^2 cm-3) 


iVn.rcs 

(10^0 cm-2) 


^out 

(10" cm) 


J0420 


45.0 


9/0.345 


22.5 


2.0 


3.0 


<1.2 


J0720 


88.4 


9/0.36 


35.4 


1.6 


2.1 


<2.2 


J0806 


87.2 


3.5/0.25 


43.6 


1.6 


2.6 


< 1.2 


J1308 


102 


5.5/0.5 


51 


1.4 


2.7 


< 1.8 


J1605 


92.6 


7.5/0.39 


23.2 


1.8 


2.0 


1.2 


J1856 


63.5 


10/0.16 


25.4 


1.8 


1.7 


< 1.5 


J2143 


104 


4.5/0.43 


10.4 


1.5 


4.6 


2.2 



4 DISCUSSIONS AND CONCLUSIONS 

In the discussion section of Kaplan et al. (201 1), they also mentioned the RCS process, which they refer 
to the paper of Lyutikov & Gavriil (2006). However, as pointed in Tong et al. (2010), the one dimensional 
treatment of Lyutikov & Gavriil (2006) can only result in up scattering. While the optical/UV excess 
of XDINSs is of "soft excess" in nature, which requires down scattering process. This is the reason 
that Tong et al. (2010) modeled the RCS process three dimensionally and employed the Kompaneets 
equation method. 

The RCS model can explain the spectral energy distributions of XDINSs in a wide parameter space. 
Meanwhile, since the surface X-ray photons are scattered by a shell of electrons, this will result in a low 
pulsation amplitude naturally. The X-ray spectra are Planck-like with no high energy tails. This may 
imply that XDINSs are possibly quark stars (Xu 2002,2009). At the same time, the possible features 
may be proton cyclotron lines or electron cyclotron lines far in the magnetosphere (Turolla 2009). 

From figure 3 in Kaplan et al. (2011) we see that for RX J2143.0H-0654 we have only two data 
points at present. Future more optical and UV data will tell us whether its spectrum is really so flat or 
not. Moreover, with more data points we can know whether there is curvature in the optical/UV spectrum 
and whether there are more spectra components in addition to the power law component. The electron 
thermal bresstrahlung component has a flat spectrum in the IR/optical range, with a flux about O.lyitJy 
for RX J2143.0H-0654. This flux is much lower than the i7-band flux upper limit (1.6/iJy at 1.65 i-im, 
Lo Curto et al. 2007). Therefore in order to verify whether the flat spectrum of RXJ2 143. Oh-0654 is due 
to thermal bremsstrahlung process, deeper optical/UV and IR observations are required. 

The electron system cools down via the bremsstrahlung process (possibly also via free-bound transi- 
tions, i.e., recombination process). Therefore, the central star has to continually heat the electron system 
and ionize it. This is similar to the "Stromgren sphere" around early-type stars (Dyson & Williams 
1980). The size of the ionized region is (3S'/47r7V2^2)^^^ (Dyson & WilHamas 1980), where S is the 
ionizing photon number flux of the central star, N^ is the electron number density (assuming fully ion- 
ized hydrogen plasma), /32 is the recombination coefficient. Similar calculations show that the size of 
the ionized region is ^ 10^° cm, to order of magnitude consistent with the outer radius used in table 1 . 
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Appendix A: FORMATION OF THE ELECTRON SYSTEM 

The electron system may be formed similarly to the "injection through ionization" process (Luo 
& Melrose 2007). Assuming a neutron star whose mass is 1.4 M0, radius 10 km and spatial ve- 



6 Tong, Xu & Song 

locity 100 km s^^, it will accrete interstellar medium or circumpulsar material at a rate M = 
2t:{GM)'^V~^Pism ^ 2 x 10® gs"^ (Frank et al. 2002), where the density of interstellar medium 
is assumed to be 10~^^ g cm^'^. Since the accretion rate is rather low, the corresponding Alfvcn radius 
Tm ~ 4 X 10^^ cm is larger than the light cylinder radius and corotation radius. Only a small fraction 
of the accreted matter can fall onto the neutron star through diffusion process, similar to the propeller 
case (Ertan et al. 2007). Since we are considering diffusion process along magnetic field lines, the dif- 
fusion properties are not affected by the presence of magnetic field. Assuming the plasma deflection 
length as the typical mean free path (Frank et al. 2002), the particle flux of the diffusion process is 
F = -D|VnisM| ^ I Aiirms "r^" ^ 10*'cm~^s~^, where D is diffusion coefficient, A the mean free 
path ~ 10^^ cm, fims the typical particle velocity ~ 10kms~^, and riisM the interstellar medium 
number density ~ 1 cm"'^. With such a diffusion flux, the corresponding accretion rate is about two 
hundred times smaller than the total Bondi-Hoyle accretion rate given above. The time scale for the 
formation of the electron system is about 4 x lO'^ yr <C Tc , where Tc is the pulsar characteristic age 
« 4 X lO^yr for RXJ2143. 0+0654. When the diffused plasma comes to the vicinity of the neutron star, 
they will diffusion across magnetic field lines and enter into the pulsar closed field regions. The plasma 
system of the RCS model may be formed. The maximum particle number density that can be confined 
by the magnetic field at Tout is (i.e., the kinematic energy density and magnetic one are in equilibrium): 
-B(rout)^/(87rfcro) ^ 10^^ cm""^ ^ Nc- The plasma system can be confined by the magnetic field. 

References 

Dyson, J. E., Williams, D. A. 1980, Physics of the interstellar medium, Manchester University Press, 

Manchester 

Ertan, U., Erkut, M. H., Eksi, K. Y, Alpar, M. A. 2007, ApJ, 657, 441 

Frank, J., King, A., Raine, D. 2002, Accretion power in astrophysics, Cambridge University Press, 

Cambridge 

Kaplan, D. L., van Kerkwijk, M. H. 2009, ApJ, 705, 798 

Kaplan, D. L., Kamble, A., van Kerkwijk, M. H., Ho, W. C. G. 201 1, ApJ, 736, 1 17 

Lo Curto, G., Mignani, R. P, Perna, R., Israel, G. L. 2007, A&A, 473, 539 

Luo, Q., Melrose, D. 2007, MNRAS, 378, 1481 

Lyutikov, M., Gavriil, F P 2006, MNRAS, 368, 690 

Ruderman, M. 2003, arXiv: astroph/03 10777 

Rybicki, G. B., Lightman, A. P. 1979, Radiative process in astrophysics, John Wiley & Sons, New 

York 

Tong, H., Peng, Q. H. 2007, ChJAA, 7, 809 

Tong, H., Peng, Q. H., Bai, H. 2008, ChJAA, 8, 269 

Tong, H., Xu, R. X., Peng, Q. H., Song, L. M. 2010, RAA, 10, 553 

Turolla, R. 2009, in Neutron stars and pulsars. Astrophysics and Space Science Library, 357, 141 

van Kerkwijk, M. H., Kaplan, D. L. 2007, ApSS, 308, 191 

Xu, R. X. 2002, ApJ, 570, L65 

Xu, R. X. 2009, J.Phys.G: Nucl. Part. Phys., 36, 064010 



